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We announce the discovery of the transiting planet CoRoT-13b. Ground based follow-up in CFHT and IAC80 confirmed CoRoT's observations. 
The mass of the planet was measured with the HARPS spectrograph and the properties of the host star were obtained analyzing HIRES spectra 
from the Keck telescope. It is a hot Jupiter-like planet with an orbital period of 4.04 days, 1.3 Jupiter masses, 0.9 Jupiter radii, and a density of 
2.34 g cm"^ It orbits a GOV star with Teff = 5945K,M, = \.Q9Mq,R, = l.Oli^o. solar metallicity, a lithium content of -I- 1.45 dex, and an estimated 
age between 0.12 and 3.15 Gyr. The lithium abundance of the star is consistent with its effective temperature, activity level, and age range derived 
from the stellar analysis. The density of the planet is extreme for its mass. It implies the existence of an amount of heavy elements with a mass 
between about 140 and 300Me. 

Key words, stars: planetary systems - techniques: photometry - techniques: radial velocities - techniques: spectroscopic 
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1. Introduction 

Transiting planets are fundamental objects for the understand- 
ing of planetary formation and evolution. Their particular ge- 
ometrical orientation allows the simultaneous measurement of 
their mass and radius, permitting a first order study of their in- 
ternal structure. Moreover, through the careful analysis of their 
passages in front of (primary transit) and behind (occultation of 
the planet or secondary transit) their host star, a characterization 
of the composition and the temperature structure of their atmo- 
sphere can be caiTied out. 



* The CoRoT space mission, launched on December 27th 2006, 
has been developed and is operated by CNES, with the contribution 
of Austria, Belgium, Brazil , ESA (RSSD and Science Programme), 
Gennany and Spain. Part of the observations were obtained at the 
Canada-France-Hawaii Telescope (CFHT) which is operated by the 
National Research Council of Canada, the Institut National des Sciences 
de rUnivers of the Centre National de la Recherche Scientifique of 
France, and the University of Hawaii. Based on observations made 
with HARPS spectrograph on the 3.6-m European Organisation for 
Astronomical Research in the Southern Hemisphere telescope at La 
Silla Observatory, Chile (ESO program 184.C-0639). Based on ob- 
servations made with the IAC80 telescope operated on the island of 
Tenerife by the Instituto de Astrofisica de Canarias in the Spanish 
Observatorio del Teide. Part of the data presented herein were obtained 
at the W.M. Keck Observatory, which is operated as a scientific part- 
nership among the California Institute of Technology, the University of 
California and the National Aeronautics and Space Administration. The 
Observatory was made possible by the generous financial support of the 
W.M. Keck Foundation. 



CoRoT is a space telescope dedicated to the study of as- 
teroseismology and the discovery of extrasolar planets by the 
method of transits ( BagUn et al.|2006) . The Coi?oT survey has 
previously detected, among other planets, the faint signal of a 



small rocky planet (CoRoT-7b; Leger, Rouan, Schneider et al 



2009) IQueloz et al.|2009[ [Bruntt et al.|2010b| ) and a giant planet 
with a temperate surface temperature (CoRoT-9b; Deeg et al. 
2010] ). 

This paper presents the discovery of the planet CoRoT-13b, 
a dense giant planet in a close orbit around a main sequence star. 
Sectionl2]describes the observations performed with the satellite. 
Ground based observations including spectroscopic characteri- 
zation of the star and radial velocity measurements are repotted 
in section [3] The derivation of the planetary parameters is de- 
scribed in section [4] A discussion of the results is presented in 
section |5] 



2. CoRoT observations 

The observations of the CoRoT field LRa02 started on 
November 16th 2008; the first transits of CoRoT-13b were dis- 
covered by the Alarm Mode and the target was oversampled (the 



standard 5 12s sampling rate was changed to 32s; see Surace et al 



2008| ) since December 9th 2008 and until the end of the obser- 
vations in March 11th 2009, gathering in total about 250000 
photometric measurements. 

The coordinates, identification labels and magnitudes of 
CoRoT-13 are given in Table [T] It is a relatively faint object for 
CoRoT (V - 15.039) and it was therefore assigned a monochro- 
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light curve of CoRoT 13b 
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Fig. 1. Raw light curve of CoRoT-13b. The target was oversampled since December 9th 2008 (HJD 2454 810.20); for visualization 
purposes we superimpose the LC binned at a sampling rate of 512s to the 32s sampled region. The jumps in the data are from hot 
pixels as described in the text. 



matic aperture mask (for a complete description of the observ- 
ing modes of the satellite, please refer to|Boisnard & Auvergne| 



2006} [Barge et al.|2008[[Auvergne et al.|2009[ ). Figure |1| shows 
the raw Ught curve (LC) of the target. One of the most signifi- 
cant environmental effects in CojRoTLCs are hot pixels: proton 
impacts that produce permanent damage to the CCD lattice (see 
[Drummond et al.|2008|[Pinheiro da Silva et al.|2008] l. On aver- 
age, 0.3% of the pixels are affected by one of these events with 
an impact of more than 1000 electrons in flux during a long run 
(Auvergne et al. 2009). The mask of CoRoT-13b has only 69 
pixels, but unfortunately this particular LC shows at least 3 im- 
pacts over the 1000 electrons threshold plus some other impacts 
of less importance. These impacts do not affect significantly the 
characterization of the target as long as the regions affected are 
properly treated ( Pinheiro da Silva et al.|2008 1. 



3. Ground-based observations 

3.1. Photometric measurements 

A complete description of the photometric follow-up of 



Coi^oT targets is given in Deeg et al. (20091. For the particular 



case of CoRoT-13b, photometric measurements of the host star 
were carried out at the 3.6m Canada-France-Hawaii-Telescope 
(CFHT) in late November 2009 confirming that the transit was 
on-target, although with a depth slightly bigger than expected. 



This is due to the contribution of the contaminants to the flux 
in the CoRoT mask. Later on, in early January 2010, new on- 
off observations were carried out at the IAC80 telescope, in 
Tenerife, confirming the previous results. 

Figure |2] shows the region of the sky around the target taken 
with the Wide Field Camera of the Isaac Newton Telescope 
(INT/WFC), in Roque de los Muchachos. The closest contam- 
inants to the main target, according to the ExoDat database 
( Deleuil et al. 12009 ), are shown in Table l2] In particular, the ob- 
ject 110839426 is completely included in the Coi?oT mask, di- 
luting the transit signal of CoRoT-13b. An analysis of the point 
spread function (PSF) of Coi^oT reveals that the main target, the 
planet-hosting star, is responsible for 89 ± 1% of the flux within 
the mask. 



3.2. Spectroscopic measurements 

The star was observed with HIRES at the Keck on Dec. 5th 
2009 as part of the NASA's key science programme in support 
of the CoRoT mission. We obtained one spectrum without the 
iodine cell and with an exposure time of 1200 sec. The spec- 
tral resolution is ~ 45 000. To determine the atmospheric param- 
eters of CoRoT-13 we analyzed the spectrum using the VWA 
software (Br untt et al.||2004| |2010a| l. We selected Unes in the 
range 5050-78 10 A. A small section of the spectrum is shown 
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Table 1. IDs, coordinates and magnitudes of CoRoT-13. 



CoRoT window ID 
CoRoT ID 
UCAC3 
USNO-A2 ID 
USNO-Bl ID 
2MASS ID 

Coordinates 



RA (J2000) 
Dec (J2000) 

Magnitudes 



Filter 

y 

v 
s" 

K* 



LRa02_E2_2165 

110839339 

170-048045 

0825-03324928 

0849-0108714 

06505307-0505112 



6h 50m 53.07s 
-5° 5' 11.26" 



Mag 


Error 


15.777 


0.077 


15.039 


0.041 


14.738 


0.027 


14.304 


0.033 


13.710 


0.021 


13.406 


0.027 


13.376 


0.038 



Notes. <"> Provided by ExoDat ( JDeleuil et al.|2009| ; "'> from 2MASS 
catalog. 



Table 2. Closest contaminants to CoRoT-13b with their respec- 
tive magnitudes and relative distances to the target. 



CoRoT ID 


distance 
(arc sec) 


B 


V 


R 


110839426 


6.58 


19.354 


18.136 


17.625 


110839769 


20.14 


15.424 


14.689 


14.395 


110838938 


21.58 


19.691 


18.412 


17.838 


110838780 


22.79 


16.927 


16.066 


15.701 


110838726 


28.34 


19.911 


18.489 


17.875 


110839832 


28.44 


17.611 


16.749 


16.388 



Table 3. Abundances of 10 elements in CoRoT-13 relative to the 
Sun. The element name, abundance relative to the Sun, and the 
number of spectral lines used are given. 



El. 



Abund. 



N El. 



Abund. 



N 



Lit 


+ 1.45 




1 


Til 


-0.06 


2 


Oi 


+0.02 




2 


Tin 


+0.05 ± 0.04 


2 


Nai 


-0.07 




2 


Cri 


-0.03 ±0.16 


4 


Mgi 


-0.07 




1 


Fei 


+0.03 ± 0.04 


103 


Si I 


+0.04 ± 


0.04 


13 


Fe II 


+0.04 ± 0.05 


11 


Cai 


+0.03 ± 


0.05 


6 


Nil 


-0.04 ± 0.05 


25 



1 arcmin 
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Fig. 4. Mean abundances for 10 elements in CoRoT-13 Keck 
spectrum. White circles correspond to neutral lines and red 
boxes stand for singly ionized lines. The yellow area represents 
the mean metallicity within one sigma error bar. 



Fig. 2. Region of the sky around CoRoT-13. Left: Image taken 
with INT/WFC with a resolution of 0.3". The position of the 
target is indicated with a white cross. Right: Image taken by 
CoRoT with the same scale and orientation; the resolution is 
around 2.32"/pixel. The area inside the black line is the mask 
used to compute the photometry; the position of the nearest con- 
taminants from the ExoDat database are indicated with different 
types of crosses. 



in Fig. [3] where the typical signal-to-noise (SN) ratio is 55. To 
determine the atmospheric parameters we use the neutral and 
ionized Fe lines and also the wide Ca lines at 6122 and 6162 A 



(see Bruntt et al. 2010a for a description). From this analysis we 
determine the parameters T^g = 5945±90K, log^ = 4.30±0.10, 
[M/H] = +0.01 + 0.07 (mean of Si, Fe and Ni), and projected 
rotational velocity v sin ; =4+1 km/s. The r^ff has been adjusted 
by -40 K based on the comparison of 10 stars with Teff deter- 
mined from both interferometric and spectroscopic methods as 
described by Bruntt et al. (j2010ai. The abundances of 10 ele- 
ments are given in Table |3] and shown in Fig. |4] The horizontal 



yellow bar is the mean metallicity. Although the metallicity of 
this star is essentially solar, it is notable that the abundance of Li 
I is +1.45 dex, in similarity to the case of another planet hosting 
star, CoRoT-6 ( Fridlund et al. 2010 ). 

The mass, radius and age of the star are calculated using 
the Starevol evolutionary tracks ( Siess| 2006 Palacios, private 
communication). The input values are the spectroscopic param- 
eters derived above (Tgff and [Fe/H]) and the proxy for the stel- 
lar density (MJ /7?*) obtained from the modeling of the tran- 
sit light-curve. The results are M, = 1.09 + 0.02Mo; /?, = 
1.01 + O.Q3Rq; which gives a corresponding surface gravity of 
logg - 4.46 ± 0.05; in perfect agreement with the spectroscopic 
value. We obtain an interval for the star's age between 0.12 and 
3.15 Gyr. 

Finally, using simultaneously the seven BVr'i'JHKs broad- 
band magnitudes reported in Table [T] and following the method 
described in iGandolfi et al.| ( |2008| l, we derived an interstellar 
extinction Ay = 0.20 + 0.10 mag and a distance to the star 
d - 1060 + 100 pc. Consistent results were obtained using the 
absolute magnitude (Allen_1973) and the intrinsi c near-infrared 
colours ( jGonzalez Hernandez & Bonifacio|2009[ ) tabulated for a 
GO V type star. 
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Fig. 3. Detail of the Keck spectrum of CoRoT-13 around the pressure-sensitive Ca n line at 6162 A. The synthetic spectrum is 
shown with a dashed line. Spectral lines used in the abundance analysis are marked with solid vertical lines while the other spectral 
lines are marked with dotted Unes. 



3.3. Radial Velocity measurements 

Precise radial velocity measurements of CoRoT-13 were ob- 
tained with the HARPS spectrograph between the nights of 
November 22nd, 2009 and February 15th, 2010 (ESO program 
184.C-0639). HARPS is a cross-dispersed echelle spectrograph 
fiber- fed from the Cassegrain focus of the 3.6 m telescope at La 
Silla Observatory, Chile ( .Mayor et al. 2003 i. Fifteen spectra with 
a spectral resolution /J ss 1 15 000 were obtained using exposure 
times of 3600 s, and setting one of the two available fibers on 
the sky in order to monitor the presence of moonlight and to ob- 
tain an optimal sky background subtraction, which is important 
for faint targets such as this. The signal-to-noise ratios per pixel 
at 5500 A of these observations range from 7.1 to 11.7. Th-Ar 
calibrations were obtained at the beginning of each night, which 
has been shown to be enough to obtain the required precision, 
due to the high stability of the instrument. 

The spectra were reduced and extracted using the HARPS 
pipeline, and the radial velocity was measured on each extracted 
spectrum by means of a weighted cross-correlation (see Baranne 
et al.|199 6) with a numerical mask corresponding to a G2 star. 
The resulting cross-correlation functions (CCFs) were fit by 
Gaussians to get the radial velocities. The measured values are 
listed in Table |4] and shown in Fig. E\ together with the best fit 
orbital solution (see below). During some of the observations, 
the star fiber was contaminated by moonlight. In those cases, if 
the peak of the CCF produced by moonlight was expected to be 
close to the measured speed of the target, a coiTection was ap- 
plied using the fiber which recorded the sky (see Bonomo et al. 
2010, in preparation). Those points are shown as white circles 
in Fig. [sjand we added quadratically 30 m s' to the uncertainty 
estimated from the CCF, in order to account for possible system- 
atic eiTors introduced by the moonlight coiTection. 

The orbital solution was found by x^ minimization, with the 
period and epoch of inferior conjunction (when radial velocity 
is zero after removal of the systemic velocity) being fixed to the 
values provided by the CoRoT ephemeris (which are calculated 
by fitting a linear regression to the center position of the indi- 
vidual transits). The eccentricity of the orbit was a free parame- 



ter at first, but since the best fit solution was compatible with a 
circular orbit at the two-cr level (the three-cr upper limit to the 
eccentricity is 0.145), we decided to fix it to e = for the de- 
termination of the rest of the parameters and their uncertainties. 
FigurelSlshows the RV measurements, phased to the CoRoT pe- 
riod, together with the best fit circular model and the residuals; 
the obtained parameters are listed in Tablel6] The resulting value 
of;t'^ is 8.2 for 13 degrees of freedom, and the rms of the resid- 
uals is 20.2 m s"', which is compatible with what should be ex- 
pected based on the median of the RV uncertainties, 21.2 ms"'. 
These facts suggest that the circular model - with the obtained 
parameters- adequately describes the available data. 

With fixed ephemeris and eccentricity set to zero, the fitting 
problem is reduced to a linear least-square minimization with 
two free parameters. The uncertainties reported in Table |6] are 
therefore estimated by means of the covariance matrix, which 
has a covariance term of 6.48 m^ s"^. However, stellar activity 
and other long-term phenomena can produce correlated noise in 
the observations and hence render the above estimation of the 
uncertainties invalid. In order to explore this we used the Prayer 
Bead method (see, for example, | Desert et al.|2009 Winn et al 



2009 1, i.e. we performed a cyclic permutation of the residuals 



of the best fit curve and fit the model again. We repeated this 
for every possible shift and measured the standard deviation of 
the obtained parameters. We also performed a similar analysis 
but randomly re-ordering the residuals rather than shifting them. 
In this way, we constructed 10000 synthetic data sets that were 
used to fit our model again. In both cases, the obtained dispersion 
of the parameters were smaller than the error bars reported in 
Tabled 

The bisector analysis for these data is shown in Figl6] where 
the uncertainty in the bisector span velocity has been set to twice 
that of the corresponding radial velocity. The bisector span ve- 
locities do not show any clear dependence with radial velocity 
values and the Pearson correlation coefficient between these two 
magnitudes is around 0.15, which is a sign of lack of coiTe- 
lation. This fact clearly indicates that the measured RV varia- 
tions do not originate from changes in the shape in the CCF 
as would be the case if the system consisted of a background 
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Table 4. Radial velocities measurements and bisector span ve- 
locities measured with HARPS. 



BJD 


RV 


ctrv 


Bis 


Moon Correction" 


-2455000 


[km/s] 


[km/s] 


[km/s] 


[km/s] 


158.7187 


22.293 


0.029 


0.027 




160.725 


22.598 


0.026 


0.038 




161.7209 


22.522 


0.053 


0.144 


0.101 


163.8279 


22.379 


0.044 


-0.019 


0.055 


166.7245 


22.311 


0.037 


-0.091 


0.219 


168.7634 


22.641 


0.042 


-0.012 


0.503 


219.7315 


22.296 


0.02 


0.036 




220.7327 


22.464 


0.022 


-0.014 




225.6511 


22.606 


0.021 


-0.01 




226.7338 


22.45 


0.021 


0.059 




228.6545 


22.452 


0.021 


-0 




237.5841 


22.624 


0.02 


0.006 




238.5617 


22.478 


0.018 


-0.014 




239.7142 


22.304 


0.019 


0.024 




243.6139 


22.311 


0.016 


-0.003 





*"' Difference between moonlight-corrected and uncorrected radial ve- 
locities. 
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Fig. 5. Upper panel: Radial Velocity measurements phased to 
the orbital period measured by CoRoT. The solid curve repre- 
sents the best fit solution. Lower panel: Residuals to the fit. In 
both panels, the white symbols indicate the measurements that 
have been corrected for moonlight contamination. 



eclipsing binary whose light were diluted by the CoRoT main 
target. Additionally, no significant changes in the measured RV 
are observed when different stellar masks are used for the corre- 
lation, which further excludes the background eclipsing binary 
scenario. 



4. Planetary parameters 

The transit parameters were determined in the same manner as 
described in the case of CoRoT-6b (Fridl und et al.||20T0| l. The 
pre-processed light curve was divided by its median, and a new 
light curve was constructed by convolving it with a fourth order 
Savitzky-Golay filter ( Press et al.|2002| l. The standard deviation 
of the differences between the measured and the convolved light 
curves is calculated and a 5cr clip is applied to remove spurious 
outliers. This clipping is iterated until no more outliers are found. 




22.30 22.35 22.40 22.45 22.50 22.55 22.60 22.65 22.70 

Radial velocity [km/s] 



Fig. 6. Bisector analysis of the HARPS data. The white symbols 
represent measurements that have been corrected for moonlight 
contamination. 



After these steps, the Ught curve is folded and a bin average is 
applied forming 5000 bins. 

We used the Mandel & AgolJ ( |2002| l model to perform the 
fit to the final, phase-folded light curve. We assumed a circular 
orbit and an absolutely dark planet at the CoRoT wavelength. 

The six adjusted parameters were: the a/R^, k - Rp\/R~i ra- 
tios, the impact parameter b - a * cos ijR^ (i is the inclination 
of the orbit), the two limb darkening coefficients m+, m_ and the 
contaminant factor, which was varied between 10-12% (corre- 
sponding to the measured value of 11 ± 1%) to allow a bet- 
ter estimation of errors. The quadratic limb darkening law is 
/(yu)//(l) = 1 - yu -H M„/v + Mfo(l - pf", and we fit the coefficients 
u+ = Ua + Uh and u^ - Ug - Ub- The best fit was found by the 
Harmony Search algorithm ( Geem et al.||2001 1, a genetic-type 
algorithm. The Icr errors were obtained from the width of the 
parameter distribution to be between y^ . and y^ ■ + I. 

A /i min ^ mm 

We calculated different solutions for the light curve of 
CoRoT-13b using different approaches in the treatment of the 
limb darkening parameters in order to understand the constraints 
that this effect puts in the determination of the final values. The 
results are presented in Table l5] In model A we left both u+ 
and M_ to be free parameters. In models B, C and D we fixed 
u+ at 0.81 (the value what we found in model A), at 0.88 (i.e. 
0.81 + Icr) and 0.74 (0.81 - Icr; where cr means the error bar 
of the coefficient found in model A) respectively. In model E we 
fixed the limb darkening to the predicted values for the star ( |Sing| 
2010). The coiTesponding ;if^-values of these solutions were not 
too different from each other, showing that limb darkening is a 
second order effect in the transit shape. The different assump- 
tions for the limb-darkening yielded quite consistent light curve 
solutions: the a/R^, k ratios are in good agreement well within 
their error bars in these five models. When we fixed u+, the un- 
certainties of the impact parameter increased by a factor of two, 
however, the eiTor bars of impact parameters found in the B, C, 
D and E solutions are overlapping with the one found in model 
A. The precise determination of the impact parameter is more 
complicated when the transit is nearly central. We learn from 
this experiment that fitting both limb darkening coefficients (or, 
what we actually did, fitting both of their combinations) yielded 
a bit more precise values. 

Since m+ = 0.81 ± 0.07 (Sing 201()]gives 0.662 + 0.022 for 
this temperature, log g and metallicity) and m_ = -0.09 ± 0.09 
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Table 6. Planet and star parameters. 



0.0030 
0.0000 




Fig. 7. Folded transit of CoRoT-13b, best fit and its residuals 
using the values of solution A in Table l5] 



( Sing|2010 gives 0.156 + 0.022) the agreement between theoret- 
ical predictions and measurements is withing 2 sigma error bars, 
which is satisfactory. We accept solution A as the definitive one 
(the fit is shown in Fig.lTli. 

5. Discussion 



5. 1 . Stellar properties 

The spectroscopic analysis of CoRoT-13 reveals a GOV star with 
an age between 0.12 and 3.15 Gyr, solar metallicity ([M/H] - 
+0.01 ± 0.07), a high relative abundance of lithium (+1.45 dex), 
and a low activity level according to the analysis of activity in- 
dicators such as the H-K Ca ii lines (where no emission is de- 
tected). 

The rate of lithium depletion of solar like stars is related to 
the age of the star and to the depth of the convective zone, as it 
is destroyed at a temperature of approximately ~ 2.5 ■ 10^ K in 
the radiative region of a star ( |Chaboyer|1998 l. Given the spec- 
tral type of CoRoT-13, we expect a lower Uthium depletion rate 
than in solar analogs ( Castro et al.|2009 1. Using the value of the 
Li I abundance, we compute a logn(LO = 2.55. From Fig. 7 of 
ISestito & Randich|2005| and with the value of the effective tem- 
perature (Tfiff - 5945 ± 90 K) we estimate the age of the star 
in the range 300 Myr to 1 Gyr, consistent with the range from 
the evolutionary models. In a recent paper, Israelian et al. ( |2009| l 
claim a lithium depletion in solar like stars with orbiting planets, 
although it is not clear that previous observations support this 
conclusion (" Melendez et al.p009b| l. CoRoT-13 is not depleted 
in lithium, albeit we call the attention to the fact that the effec- 
tive temperature of this star is slightly higher than the upper limit 
for depletion given in Israelian et al.^(2009| . 

The V sin / value indicates a ro tational period of the star of 
around 13 day^ Gyrochronology ( Barnes|200 7') could be used 
as an age estimator. Using the improved gyrochronology rela- 
tions from Mamajek & Hillenbrand ( |2008 1 we derived a gy- 
rochronologic age of 1 .66 Gyr well within the range of age given 
by evolutionary models. No emission feature is seen at the bot- 
tom of the Ca ii H and K lines nor in the Ha line showing that the 
star belongs to the inactive population with log R'^K < -5.0. 
We thus didn't derive any chromospheric age for the star. Other 



Ephemeris 


Planet orbital period P [days] 


4.035 190 ±0.000030 


Primary transit epoch T,, [HJD-2 450 000] 


4 790.8091 ±0.000 6 


Primary transit duration dt^ [h] 


3.14 ±0.01 


Results from radial velocity observations 




Orbital eccentricity e 


(fixed) 


Radial velocity semi-amplitude A^ [ms"'] 


157.8 ±7.7 


Systemic velocity Vr [kms"'] 


22.4536 ± 0.0060 


0-C residuals [ms-'l 


20.2 


Fitted transit parameters 




Scaled semi-major axis ajR, 


10.81 ±0.32 


Radius ratio k = Rp/R, 


0.0909 ±0.00f 4 


Quadratic limb darkening coefficients" u+ 


0.8f ±0.07 


K_ 


-0.09 ± 0.09 


Impact parameter* b 


0.374 ± 0.054 


Deduced transit parameters 




M'J'/R, [solar units] 


1.0f4± 0.030 


Stellar density p, [ g cm"^] 


1.468 ±0.131 


Inclination / [deg] 


88.02!°i^ 


Spectroscopic parameters 




Effective temperature T^g [K] 


5 945 ± 90 


Surface gravity log g [dex] 


4.30 ±0.10 


Metallicity [Fe/H] [dex] 


O.Of ±0.07 


Stellar rotational velocity v sin ; [ km s" ' ] 


4± 1 


Spectral type 


GOV 


Stellar and planetary physical parameters f 


•■om combined analysis 


Star mass [Mg] 


1.09 ±0.02 


Star radius [Rq] 


1.01 ±0.03 


Surface gravity log g [dex] 


4.46 ± 0.05 


Age of the star / [Gyr] 


0.12-3.15 


Distance of the system [pc] 


1 060 ± 100 


interstellar extinction Ay [mag] 


0.20 ±0.10 


Stellar rotation period Pjot [days] 


13- 


Orbital semi-major axis a [AU] 


0.05 10 ±0.0031 


Planet mass Mp [Mj^p]^ 


1.308 ±0.066 


Planet radius R^ [R]„p]^ 


0.885 ±0.014 


Planet density pp [gcm^^^] 


2.34 ±0.23 


Planet surface gravity log g [dex] 


3.62 ± 0.03 


Average surface temperature'' Tp [K] 


~ 1700 



' a lower limit, as the value of sin / for the spin axis of the star is 
unknown. 



Notes. <"' /(^)//(l) = I - p + Uap + util - pf, where /(I) is the 
specific intensity at the center of the disk and p = cos y, y being the 
angle between the surface normal and the line of sight; u+ = Ua + iib 
and M_ = Ua — U),. *' b = ^^^^ ''' Radius and mass of Jupiter taken 
as 71 492 km and 1.8992 x lO^" g, respectively (Lang 1999|l. <'" Zero 
albedo equilibrium temperature for an isotropic planetary emission. 



GOV stars are found with similar rotation rates and low activity 
levels ( |Noyesetal.|I984| i. 

We have looked for signs of stellar rotation in the light 
curve (LC) observed by CoRoT to make a comparison with the 
clear signs of spot modulation found in the cases of CoRoT-2b 
(jLanza et al.f2009b), CoRoT-4b (jLanza et al.||" 2009a), CoRoT - 
6b (Fridlund et al. 2010) or CoRoT-7b (Lan za et aL||20T0l l. 
The Lomb-Scargle periodogram of the LC, once the planetary 
transits have been removed and the hot-pixel events have been 
treated, shows indeed a significant broad peak around 77 days; 
but not any significant peak at the expected rotational frequen- 
cies. The 77 days period is comparable with the length of the 
run (115 days), so it might be that we are observing an irregular 



J. Cabrera et al.: Transiting exoplanets from the CoRoT space mission 

Table 5. Values of the adjusted parameters in the modehng of the transit of CoRoT- 1 3b in the different approximations described in 
the text. 



Parameter 


A 


B 


C 


D 


E 


ajR. 


10.81 + 0.32 


11.22 ±0.42 


11.15 ±0.40 


10.98 ± 0.48 


10.89 ± 0.40 


k = Rp/R. 


0.0909 ±0.0014 


0.0900 ±0.0016 


0.0896 ±0.0014 


0.0912 ±0.0012 


0.0914 ±0.0011 


b - a cos ijR, 


0.374 ± 0.054 


0.264 ± 0.097 


0.271 ±0.092 


0.349 ±0.091 


0.385 ± 0.070 


inclination (deg) 


88.01!«i 


88.65!°:^^ 


88.61!«: 


88.or°5^ 


87 97+0.43 


u+ 


0.81 ±0.07 


0.81 (fixed) 


0.88 (fixed) 


0.74 (fixed) 


0.662 (fixed) 


U- 


-0.09 ± 0.09 


-0.09 (fixed) 


-0.24 ± 0.32 


-0.01 ±0.25 


0.156 (fixed) 


third light 


0.11 ±0.01 


0.11 ±0.01 


0.11 ±0.01 


0.11 ±0.01 


0.11 ±0.01 


^2(«) 


1.00000 


1.00039 


0.99945 


1.00077 


1.00130 



Notes. '"' Normalized to the value of solution A 



pattern that, in the observing window of 1 15 days, has a typical 
timescale of variation of 77 days that may mimic an harmonic 
oscillation. The only reliable information on the nature of this 
particular signal is its characteristic amplitude of 0.5% and its 
characteristic timescale of 77 days; one has to be extremely cau- 
tious when interpreting its nature. Two immediate possibilities 
are stellar activity and an instrumental residual signal. In the lat- 
ter case, stray light or other instrumental effects such as temper- 
ature fluctuations should affect a region of the CCD (if not all), 
instead of a single target. None of the targets in the neighbor- 
hood of CoRoT-13 show a similar pattern. Other environmental 
features (such as hot pixels) although present, have a completely 
different behavior both in amplitude and in timescale and are not 
likely to be responsible for the signal. The spectroscopic analy- 
sis shows that CoRoT-13 is a quiet dwarf star A GOV star might 
show spot modulation with a characteristic timescale similar to 
the rotational period of the star as well as long term variations 
in timescales of several years ( Bahunas et al.|[T995] ), but not in 
the timescales considered here. On the other hand, slow rotating 
giant stars might show a pattern of variability with similar am- 
plitudes and characteristic timescales as the ones revealed by the 
periodogram. We therefore conclude that the modulation mea- 
sured in the LC of Coi?oT- 13 is due to a background contami- 
nant and not to the main target. 



In a recent paper, Lanza (2010 1 links the presence of giant 



planets and the angular momentum loss of their host stars, de- 
tecting that giant planet hosting stars with Teff > 6QQQK tend to 
be in a n/Q a; 1,2 orbital period/rotational period synchroniza- 
tion with their respective planets. The effective temperature of 
CoRoT-13 lays right below this lower limit and indeed the rota- 
tional period obtained from the spectroscopic analysis (around 
13 days) is longer than the expected 2:1 resonance or the syn- 
chroiuzation; although the rotational period is not sufficiently 
well constrained (the relative uncertainty is around 35%) and 
the inclination of the spin-axis of the star is not known. A study 
of the Rossiter-McLaughlin effect of this planet will provide the 
relative angle of inclination of the planet and the star, which is 
an important parameter in the model proposed by Lanza (2010[l. 



5.2. Planet interior 

Figure [8] compares CoRoT- 1 3b to other transiting exoplanets 
in a mass-radius diagram. Although it is not the densest ob- 
ject known so far (both super-Earths and brown dwarfs may be 
denser), it is clearly extremely dense for its mass. This is con- 
firmed by a combined modeling of the star and planet evolu- 
tion (see|Borde et al.|2 010; Mor el & Lebreton|2008t [GuiUot &| 
|Morel|1995| l shown in Fig. [9] The small planetary radius derived 



from the transit photometry and spectroscopy can be explained 
only by advocating the presence of a truly considerable amount 
of heavy elements in the planet. When fitting the stellar effec- 
tive temperature and density within their 1-sigma error bars (see 
Table [6] and red area in Fig. l9]l, between about 140 and 300M(b 
of heavy elements are required to reproduce the measured plane- 
tary size. When fitting the stellar parameters only within 3cr (yel- 
low area), still at least lOOMjg of heavy elements are needed. In 
our planetary evolutions calculations, we assumed all heavy ele- 
ments to be grouped into a well-defined central core, surrounded 
by a solar-composition envelope, with no added sources of heat. 
The possibility that these heavy elements may be at least partly 
mixed in the envelope is not expected to change these num- 
bers significantly (Guillot_2005 , Ikoma et al._2006, Baraffe et aE] 
i2008l ). 

This extremely high amount of heavy elements is surprising 
for several reasons: first, this is probably a record - HD 149026b, 
a Saturn-mass planet, has been known to possess about 60 - 
70Me of heavy elements (Sato et al."2005; Ikoma e t al.|2006l l. 
Other objects, such as OGLE-TR-56b and OGLE-TR-132b were 
also later shown to have close to lOOMg in heavy elements 
(GuiUot et al. 2006| l, but no planetary-mass object was yet shown 
to possess more heavy elements. An exception may be HAT-P- 
2b, initially thought to have Mz - 200 to 600Me (BaraffeetaT 



2008; Leconte et al..2009j , but a revision of the stellar parame- 
ters ( Torres et al.|2008 l yields much smaller Mz values for this 
object. 

Second, planet formation models do not predict the existence 
of such dense objects (Mord asini et al.|2009| l. This is because a 
protoplanetary core growing beyond a few tens of Earth masses 
rapidly captures any surrounding hydrogen and helium, and that 
the accretion of planetesimals is suppressed by the growth of 
the protoplanet beyond Saturn's mass. Giant impacts need to be 
invoked for a further capture of a significant mass of heavy ele- 
ments ( [Ikoma et al.|2006| l. 

Third, CoRoT-13b is a counterexample outside the correla- 
tion between stellar-metallicity and planetary mass in heavy el- 
ements HGuillot et al.pOOel [Burrows et al.|2007||Guillot|2008| l. 
One interesting possibility however is related to the high lithium 
abundance of the star. If the metallicity of stars with planets in 
fact tell us about the late accretion of circumstellar gas filtered 
of their heavy elements by planet formation ( [Melendez et al. 
2009b; R amirez et al.|2009[[Nordlund|2009l ), it raises flie possi- 
bility that the CoRoT-13 system was in fact metal-rich, but that 
a heavy-elements poor, lithium-rich rare last burst of accretion 
modified the chemical properties of the star's thin outer convec- 
tive zone. 
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Fig. 8. Position of CoRoT-13b (square) among tlie other transit- 
ing planets in a mass-radius diagram. 




Fig. 9. Age (in Ga= 10^ years) versus transit radius of CoRoT- 
13b (in Jupiter units, lRjup=71 492 km). The colored area cor- 
respond to constraints derived from stellar evolution models 
matching the stellar density and effective temperature within 
a certain number of standard deviations: less than Icr (red), 
2cr (blue) or 3cr (yellow). The curves are evolution tracks for 



CoRoT-13b (assuming M = 1.308Mjup, 
various models as labelled 



T,„ = 1700 K), with 



5.3. Thermal losses 

We estimate the thermal mass loss of CoRoT-13b by using the 
method and formulae described in |Lammer et al.| ( [2009| l, and we 
find a negligible escape rate which did not influence the planets 
mass over its history. The reason of the negligible thermal mass 
loss of CoRoT-13b is the planets compactness and high density 
of about 2.34 g cm"-'. 

5.4. Occultation of the planet by the star 

In the visible wavelength range, occultations of planets behind 
their host stars have to date been detected only in a few cases: 
Co rot-lb ([Snellen et 31120091 |Alonso et at |2009a| l and CoRoT- 
2b ( |Alonso et al.|2009b|1snellen et al.|2010|l in data from CoRoT 
and HAT-P-7 in data from Kepler ( [Borucki et al.|20 09). All ai-e 
inflated giant planets {Rp > 1 .4Rjup) on very short orbits, which 
means that that they have large surfaces favoring a large flux of 
reflected light and additionally they are very hot (Tgff > 2000K), 



Fig. 10. Expected amplitude of the occultation of the planet as 
a function of the wavelength for two extreme values for the heat 



redistribution factor for CoRoT-13b (as defined by Cowan & 
Agol 201()j e = 1: uniform heat distribution over the planet; 
e = no distribution) and two values of the Bond albedo. The 
gray shaded area is the wavelength window observed by CoRoT. 
The black arrow is the upper limit of Coi?oT observations. 



light curve of CoRoT 13b 
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Fig. 11. Raw light curve of CoRoT-13b between the 
observed transits number 8 (HJD=2 454 819.06) and 9 
(HJD=2454 823.09). a) is a detail of the 8th transit, b) is a detail 
of the event at HJD=2 454 82 1 .03 which mimics the transit of a 
long period planetary companion. 



which favors thermal emission in the visible regime. Using the 
definition from Cowan & Agol (2010), CoRoT-13b, in the most 
favorable case (zero albedo, zero heat redistribution factor; up- 
per solid line in Fig. [T0|, should have an average surface tem- 



perature around 1700K. The emission of the planet in the visible 
range observed by CoRoT is therefore dominated by reflection, 
meaning that eclipse amplitudes larger than 2.5 ■ 10"^ cannot 
be expected, even with a very high Bond albedo of A^ = 0.5 
(Fig. 10 dotted lines). Such a signal is however not detectable in 



the current data, as the measured scatter in the folded and binned 
light curve is on the order of 10 '*. 

5.5. Constraints on the presence of additional planets. 

The region of the light curve (LC) shown in Fig. [TTl between 
the 8th and the 9th transits observed by CoRoT shows a feature 
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Fig. 12. Parameter space of the additional planets in the system 
which can been discarded with current HARPS measurements. 
The region above the curve is discarded by RV measurements, 
while planets located in the region below the curve would pro- 
duce an undetectable signal. 



which mimics a single planetary transit not related to the pattern 
of transits of CoRoT-13b. Almost in the middle point of the two 

regular transits we find a short feature, roughly 1 .5h long and 
1 % deep which in principle could be due to a planetary transit 
from an additional planet in the system. This feature does not re- 
peat in the LC so that a hypothetical companion would need a pe- 
riod longer than the 80 days that remained until either end of the 
observations. Nevertheless, the duration of a planetary transit is 
linked to the period of the orbit ( Seager & Mallen-Ornelas|2003| l 
and for an 80 days orbit we expect a transit around 8 hours long. 
Of course, some corrections have to be done in the case the orbit 
is eccentric ( Carter et al.||2008] |Kipping||2008| l; but in our case 
we would have to admit an eccentricity of 0.94 (and a favorable 
orientation) for the hypothetical planet to have a 1.5h long tran- 
sit in such a long period orbit. However, this hypothetical planet 
would have a periastron passage inside the orbit of CoRoT-13b. 
There are studies of the planetary three body problem with high 
eccentricities (see for example Beauge & Michtchenko 2003; 
Michtchenko et aL]|2006] l, but to the knowledge of the authors 
there is no possible justification for the dynamical stability of 
one system formed by a close-in Jupiter-like planet in a circular 
orbit and a highly eccentric Jupiter-sized companion. 

Moreover, there is no hint of the drift produced by any addi- 
tional planet in the scatter of the residuals of the circular radial 
velocity (RV) analysis. The region of the parameter space of an 
additional planet in a circular orbit which can be discarded based 



on RV data is shown in Fig 12 We note that for periods below 



85 days (the span of HARPS observations), only planets in cir- 
cular orbits with masses below ~ 0.45 Mjup are not discarded by 
RV data. The mass limit for eccentric orbits is even lower in this 
range of periods, as the amplitude K of the RV movement grows 
as(l-e2)-i/2. 

Given these facts, the most simple explanation is an environ- 
mental effect that mimics a transit. 



6. Summary 

We have reported the discovery of CoRoT-13b, a transiting giant 
planet orbiting the star CoRoT 1 10839339. 



The spectroscopic analysis reveals that CoRoT-13 is a GOV 
star with T^ff = 5 945K, M. = 1.09Mo, R, = 1.01 + 0.03Rq, 
solar metallicity ([M/H] - H-0.01 ± 0.07) and a high relative 
abundance of lithium (+1.45 dex). The evolutionary tracks con- 
strain the age of the star between 0.12 and 3.15 Gyr. The results 
of the study of the lithium abundance, the v sin / and the activity 
level (which is low both in the spectroscopic and in the photo- 
metric analysis) are consistent with the spectral type and the age 
range of the star; although the lower limit of the interval is not 
favored due to the low activity level measured. 

CoRoT-13b is the highest density jovian planet among its 
relatives with masses between 0.2 and 2.5 Jupiter masses. Its 
extreme density implies the existence of a significant amount, 
between 140 and 300M®, of heavy elements in the planet, which 
is contradiction with the expected link between Mz of a planet 
and the metalliciy of its host stars ( [Guillot et al.|2006 Burrows 
|etal.|2007t|Guillot|2008D . 

There is no hint of any other massive companion in the sys- 
tem, which is quite a common case for hot Jupiters. 

We have discussed the consequences of this characteriza- 
tion and we show evidence that this particular target provides 
important information for the knowledge of the interactions be- 
tween the structure of stars, the formation of planets, and their 



joint evolution (Ikoma et al. 2006, Amml er- von Eiff et al. 
'Melende zet al.||2009a[lNordlund||2(J09l |Ramirez et al.| 
Lanza 20To][ 



2009 



2009 
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